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Transfer-matrix methods are used to study the probability distributions of spin-spin correlation 
functions G in the two-dimensional random-field Ising model, on long strips of width L = 3 — 
15 sites, for binary field distributions at generic distance R, temperature T and field intensity 
ho- For moderately high T, and ho of the order of magnitude used in most experiments, the 
distributions are singly-peaked, though rather asymmetric. For low temperatures the single-peaked 
shape deteriorates, crossing over towards a double-^ ground-state structure. A connection is obtained 
between the probability distribution for correlation functions and the underlying distribution of 
accumulated field fiuctuations. Analytical expressions are in good agreement with numerical results 
for R/L > 1, low r, ho not too small, and near G = 1. From a finite-size ansatz at T = Tdho = 0), 
ho 0, averaged correlation functions are predicted to scale with L^ho, y = 7/8. From numerical 
data we estimate y = 0.875 ± 0.025, in excellent agreement with theory. In the same region, the 
RMS relative width W of the probability distributions varies for fixed R/L = 1 as W ^ ho f{Lho) 
with K ~ 0.45, u ~ 0.8 ; f{x) appears to saturate when x — > oo, thus implying W ^ ho in d — 2. 
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I. INTRODUCTION 

It is by now well established that the space dimension- 
ality d = 2 is the lower critical dimension of the random 
field Ising model (RFIM) in agreement with the 

early domain- wall picture of Imry and Ma Q. Thus, 
as usual for a borderline dimensionality, details of two- 
dimensional behaviour are rather intricate. The diver- 
gence of the low-temperature correlation length as the 
field intensity approaches zero is apparently anomalously 
severe ||]. This is at least partly responsible for diffi- 
culties encountered in the application of normally very 
powerful numerical techniques to this problem. In par- 
ticular, transfer- matrix (TM) methods have been used, 
either for fully finite [^-^ or semi-infinite ^ geometries. 
TM calculations have usually focused upon the structure 
factor, as obtained from suitable derivatives of the par- 
tition function. The correlation length is then derived 
from the structure factor, under assumption of specific 
scaling forms results thus far have been at least in 

qualitative agreement with theoretical predictions ||^. 

Many recent studies of the RFIM, both m d = 2 
and 3, have concentrated on zero-temperature properties, 
as an exact ground-state algorithm first applied some 
time ago jH has been revisited |jll|.|l2|. In our earlier 
work p3| , |14[ , where a domain-wall scaling picture was 
developed for bar-like systems in general d, numerical 
support for theory was provided in d = 2, T = by 
a version of such algorithm adapted to strip geometries. 
For T ^ we relied on a TM treatment of the free energy, 
again on strips. 



Here we deal directly with probability distributions of 
spin-spin correlation functions, calculated by TM meth- 
ods on semi-infinite (strip) systems. Interest in proba- 
bility distribution functions has increased recently, re- 
garding extensive quantities in critical disordered sys- 
tems. This is in line with the growing realisation that 
lack of self-averaging tends to be the rule, rather than 
the exception, e.g. for susceptibilities and magnetisa- 
tions in such systems ||l^ , implying that the width of the 
associated probability distributions is a permanent fea- 
ture that does not trivially vanish with increasing sam- 
ple size. In the present case, lack of self-averaging does 
not come as a surprise, as correlation functions are not 
extensive [Q, so the usual Brout argument is not 
expected to apply. Also, in d = 2 the random field moves 
the second-order transition to T = 0, so the d — 2 RFIM 
is off criticality at any T > 0; experimental manifesta- 
tions of microscopic features of the d — 2 RFIM come 
indirectly through (sample-averaged) non-critical prop- 
erties [18 -20 . Indeed, consideration of the crossover be- 
haviour in the vicinity of the zero-field, pure-Ising, crit- 
ical point provides interesting information, as shown in 



Section IV 



In what follows, we first discuss the ranges of spin- 
spin distance R, temperature T and random-field in- 
tensity ho for which the statistics of correlation func- 
tions display the most interesting features, and illustrate 
our choices with simple examples. We then turn to the 
connection between field- and correlation-function dis- 
tributions, and show how, in suitable limits, one can ex- 
tract the latter from the former. Next we study the line 
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T = TciJiQ = 0),ho — > 0, and use correlation functions to 
extract information on scaling behaviour corresponding 
to the destruction of long-range order by the field. A 
final section summarizes our work. 



II. NUMERICAL TECHNIQUES AND 
PARAMETER RANGES 

We calculate the spin-spin correlation function G{R) = 
(ctqCtIj), between spins on the same row (say, row 1), and 
R columns apart, of strips of a square lattice of ferro- 
magnetic Ising spins with nearest-neighbour interaction 
J = 1, of width 3 < L < 15 sites with periodic boundary 
conditions across. This is done along the lines of Sec. 
1.4 of Ref. [Q, with standard adaptations for an inho- 
mogeneous system |2^]. The strip widths used are those 
manageable on standard workstations, without unusually 
large memory or time requirements; as the main overall 
advantage of TM calculations (against e.g. those on fully 
finite, L x L systems), is that monotonic trends set in for 
relatively small strip widths |21|, the upper bound on 
L does not significantly constrain our analysis. It does, 
however, matter for the values of R used, since the inter- 
esting range of R/L is around one, where the transition 
between d = 1 and d = 2 behaviour takes place. 

Different sorts of averaging are involved in this case. 
For a given realisation of the site-dependent random 
fields, one has (for sufhciently low field intensity) a 
macroscopic ground-state degeneracy |l^. TM methods 
take into account the Boltzmann weights of all possible 
spin configurations, so they scan the whole set of avail- 
able ground states for a given realisation of quenched 
disorder. One must then promediate over many such re- 
alisations, which is done as follows. At each iteration of 
the TM from one column to the next, the random- field 
values h are drawn for each site from the binary distri- 
bution: 
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P{h) = -[5{h~ho) + 6{h + ho)] 



(1) 



By shifting the origin along the strip and accumulating 
the respective results, one can produce normalized his- 
tograms, P{G), of occurrence of G{R). With typical 
strip lengths N = 10^ columns, we generate 10"* — 10^ 
independent estimates of G{R) for R in the range 5—15 
which corresponds to R/L ^ 1, as explained above. 

In our previous study of the unfrustrated random-bond 
Ising model [p3| , the probability distribution function of 
correlations was expected to be log-normal for strictly 
one-dimensional systems [|l6|. This led us to a picture 
where, for strip width L and spin-spin distance R, the 
distribution would evolve perturbatively away from log- 
normal with increasing L/R. Thus, there we used loga- 
rithmic binning for the histograms of occurrence of G{R): 
a convenient interval of variation of lnG{R) was divided 
into, usually, 10'^ bins, each particular realization being 
assigned to the appropriate bin. As a similar starting 



point is not available here, and negative values of cor- 
relations may occur, we have resorted to a simple linear 
choice, dividing the whole [—1,1] interval of variation of 
G{R) into (again, usually 10'^) equal bins. 

The temperature and field intervals of interest are 
broadly circumscribed because spin-spin correlations are 
induced by the ferromagnetic (unit) interaction. Thus 
one must keep to values of T and ho that are not suf- 
ficient to render the coupling negligible; rough bound- 
aries, to be refined next, are Tc{ho = 0) = 2.269 . . . and 
hociT = 0) = 4 (above this latter value each spin always 
obeys the local field). 

We have found Tiow — 0.6 to be low enough to dis- 
play ground-state effects rather prominently. Recall that 
strictly atT = correlation- function histograms are triv- 
ial double-J peaks at G{R) = ±1; this reflects the frozen- 
domain structure of the ground state, which is best in- 
vestigated directly as done by others 1 1 1 1^ . Conversely, 



here we wished to investigate departures from the double- 
S shape, induced by increasing T. On the other hand, 
Thigh = 2.0 is high enough so that field fiuctuations (in 
the range of ho spelt out in the next paragraph) have 
mainly a perturbative effect. 
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FIG. 1. Normalized histograms P{G) of occurrence of G 
for strip width L = 5, length A'' — lO'' columns, R = 15, 
ho = 0.5, and T = 0.6 and 2.0. Binwidth 2 x 10"^. Vertical 
axis has linear scale in (a) and logarithmic in (b), the latter 
in order to emphasize values occurring with low frequencies. 



Experimental studies ||l^-[20|] of d = 2 dilute Ising an- 
tiferromagnets in a uniform field H (argued by Fishman 
and Aharony |24| to be equivalent to the RFIM) con- 
centrate on i?- values corresponding to /iq ^ 0.1 — 0.2 in 
Eq. (|l|) , enough to cause significant departures from zero- 
field behaviour. Higher fields /iq ^ 1 are convenient to 
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enrich domain statistics in simulations of fully finite sys- 
tems, as they reduce low-temperature domain sizes and 
increase degeneracy [|ll|,^. However, already for ho = 
0.5 the histograms of correlation functions were found to 
be utterly distorted (compared to a paradigm of single- 
peaked structures with reasonably-defined widths), so as 
to be intractable in terms of a simple description with 
few parameters. This echoes the experimental obser- 
vation for Rb2Co.7Mgo.3F4, that "...applied fields very 
much less than the Co^+ molecular fields ... have quite 
drastic effects" Fig. |l| illustrates the point. 

Nevertheless, for /iq^ 0.1 — 0.15 and high T, the over- 
all picture stays very close to that depicted in Fig. |, with 
the following main features: (i) a clearly-identifiable sin- 
gle peak, below the zero-field value Go = G{ho — 0); (ii) 
a short tail below the peak and a long one above it, such 
that (iii) all moments of order > of the distribution are 
above Go- In Figure |^ we show the zcroth (exp(lnG')) 
and first ((G)) moments. 
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FIG. 2. Normalized histogram P{G) of occurrence of G 
for strip width L = 5, length A'' = lO'' columns, R — 15, 
ho = 0.05 and T = 2.0. Binwidth 2 x 10"^. Vertical bars 
in inset located respectively at: Go (full line); exp(lnG) 
(short-dashed); (G) (long-dashed). 



This scenario breaks down for low temperatures, as 
G{ho — 0) becomes close to the upper limit of unity, for 
R/L ~ 1 and the strip widths within reach. However, 
this latter regime can be understood in terms of a di- 
rect connection between field- and correlation-function 
distributions, described in Section [II below. 



For fixed R, small ho and high T, Figure || shows the 
typical evolution of distributions against L. Note that, 
with ho — 0.15, the single-peak structure shows early 
signs of fraying. 
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FIG. 3. Normalized histograms P{G) of occurrence of G 
for strip widths L = 5, 9 and 13; length = 10^ columns, 
R = 10, ho = 0.15 and T = 2.0. Binwidth 2 x 10"^. 



One can see in Figure ^ a narrowing effect with in- 
creasing L. This is quantitatively depicted in Figure U, 
for which the use of R/L on the horizontal axis is in- 
spired in usual ideas of finite-size scaling, and has proven 
fruitful in our earlier study of random-bond systems [E3| . 
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FIG. 4. Double-logarithmic plot of RMS relative widths 
W of distributions against R/L. Strip widths L — 5, 9, 
13; R — 10, 15; ho — 0.05 (triangles) and 0.15 (squares). 
Straight lines are least-squares fits to each set of data (see 
text). T = 2.0. 
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Though the RMS relative width W = ((G^) - 
{G)^y^^/{G) appears to approach zero for small R/L 
(where, as argued in Ref. d — 2 behaviour should 
show up), we note that: (i) the data display non- 
monotonic jumps even for fixed /iq; and (ii) power-law 
fits of W against [R/LY show a rather strong depen- 
dence of X on ho (for data in Figure ^ one has x ~ 1 for 
ho = 0.05, and 2 for /iq = 0.15). These facts indicate 
that, from consideration of the above data alone, we are 
not in a position to conclude that W ^ Q as the true 
d — 2 regime is approached. Indeed, in Section IV below 



a different analysis, at fixed R/L, strongly suggests that 
the widths do not vanish in the two-dimensional limit. 



III. DISTRIBUTION OF G FROM FIELD 
DISTRIBUTION 

An important question is how the field distribution 
gives rise to the distribution for the correlation function 
G(i?) = (ctoct/j) (at specific separation R). 

A scenario worth exploring is the following: the prob- 
ability distribution P{G) for G{R) arises from a dis- 
tribution of characteristic scales ^, related to G(i?) via 
G(i?) ~ exp(— i?/^), with ^ distributed according to some 
distribution. This last probability distribution has then 
to be related to the field distribution. At low tempera- 
tures a domain picture might provide that relationship: a 
distribution of domain sizes arises from the distribution 
of fields aggregated over each domain, e.g. by minimis- 
ing energy (or free energy) along the lines of Ref. [^ , 
but generalised to consider specific field configurations, 
with their associated probability (the free energy min- 
imisation may make such an approach applicable up to 
temperatures of order Tc(h = 0) ). 

The simplest such scheme uses a common domain size 
^, over which the total field is h = xho^/^ with p{x) — 
/^/\/27r . This is the distribution of aggregated fields 
on a domain, arising from the independent distributions 
i[(5(/i — ho) + 5{h -\- ho)] of fields hi at each site i.. Then 
minimising the free energy per unit length (for the T = 
problem) gives 



e = e(-x) = 



(2) 



Hence, from the probability distribution p(x), there arises 
a probability distribution for ^(x), and via that a prob- 
abihty distribution P{G) for G(i?) - exp(-i?/C). The 
result is 



P{G,T = 0) = 



2J^L 
hlR 



1/2 



G "5 



2 L 



2^(lnl/G)i/2 



(3) 



The important parameter in this zero-temperature de- 



scription is A 



— 2J -L 



(which is of order one for i? = 15, 



The T generalisation of such pictures involves 
the entropic contribution — T(S'o -I- ^i) to the free en- 
ergy F, which includes a contribution from positioning 
of domain walls {—TSq) (see Ref. |l3j, but still allowing 
for probabilities of specific field configurations) and also 
one from the random-walk-like wandering of the domain 
walls {-TSi). 

These entropies are (using the simplest picture of a sin- 
gle C(x)) ^0 = ks ln^{x)/^{x) (using reduction valid for 
^{x) large), and Si = fcs ln/i-'"/^(a;) = kB{L/^{x))\n fi, 
with /i ~ z — 1, z — lattice coordination number. Min- 
imisation of F per unit length then gives 







hox^/LS, 



1 



(Lln^ + ln^- 1) 



(4) 



2JL JL 
The variable x is again distributed with the domain- 
aggregated field distribution p{x) = e"^ /^/V2tt which, 
via Equation (jj) then provides the distribution of ^ and 
finally the distribution for G ~ e~^^^ (along the general 
lines indicated above) . Different pairs of terms dominate 
Equation (^) in different regimes of /iq, T and L. Of 
special interest to us are the first-order low-temperature 
corrections. An approximate treatment of Equation (|^), 
valid for G near 1, gives 

P(G,T) cx P(G,T = 0)(lnl/G)-'*'^«'^/"'^ , (5) 

with P{G,T — 0) given by Equation (||). Apart from 
weakly L-dependcnt normalization factors, one should 
have 

P'(G) = P(G,T)G-"(ln(l/G))'3 = const , (6) 
where a = {2J^L/hlR) - 1, /3 ee 1/2 -I- AksT/JL . 
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FIG. 5. Histograms of occurrence of G for T — 0.6, 
ho = 0.5, R — 15, near G = 1; L — 3 (triangles), 5 (squares) 
and 7 (hexagons). Binwidth 2 x 10"'*. Full symbols: normal- 
ized histograms, P{G). Empty symbols: P'{G), Eq. (H). 
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In Figure || we check Equation ^ for T = 0.6, Hq = 
0.5, R — 15 and L = 3, 5 and 7. Use of narrow strips 
(i.e. R/L > 1) and high fields is important in order 
to produce broad distributions in the fow-temperature 
regime concerned. One sees that, indeed, the strong G- 
dependence of P{G, T) near G — 1 can be essentially 
accounted for by the factors in Equation (O). 



SCALING NEAR ZERO-FIELD CRITICAL 
POINT 



IV. 



According to theory ||5|,yj,|0||2|j25|l , the scahng be- 
haviour of the RFIM depends on the variable /iQlij""^ 
where is the random-field intensity and t ~ {T — 
Tc{ho))/Tc{ho) is a reduced temperature. For d > dc — 2 
Tc{ho) is the field-dependent temperature at which a 
sharp transition still occurs; it turns out that even in 
d = 2 the dominant terms still depend on the same 
combination, where now "Tc(/io)" denotes a pseudo- 
critical temperature marking, e.g., the location of the 
rounded specific-heat peak. This is true except for the 
d = 2 specific heat (which does not concern us directly 
here), where In /iQ-dependent terms also play an impor- 
tant role [^j26|. Further, it is predicted (2|] that the 
crossover exponent </) = 7, the pure Ising susceptibil- 
ity exponent. In d = 2, specific heat and neutron- 
scattering |20j data are in good agreement both with the 
choice of scaling variable as above, and with the exactly 
known 7 = 7/4 . 

Here we propose a direct check of scaling, as follows. 
For ho 0, near Tco = Tc{ho = 0), one expects 
"T,(/io)"= T,o-chl^^. Hence, h-"" ^ ^ hf * {T - T,^) 
apart from a small, finite shift. Setting T ^ Tcq and 
making the usual finite-size scaling ansatz t L^^/'^ 
with the pure Ising value v = 1 (this latter assumption 
to be verified), one obtains that the (finite-size) scaling 
variable at T = T^q must be 



2v 



iT = T,o,ho ^ 0) 



(7) 



with 0/2z/ = 7/8 in c? = 2. This implies that the corre- 
lation length related to the decay of ferromagnetic spin- 
spin correlations diverges along this particular line as 



y = ^/2iy. 



(8) 



From standard finite-size scaling ||2^, the correlation 
functions for distance R, strip width L, t = T — Tco — 
and random-field intensity ho are then expected to be- 
have as 



G{R, L, t = 0, ho) = L-'^ r{R/L, L^/iq) 



(9) 



In Figure ^ we show, for fixed R/L — 1, the scaling plot 
thus suggested, where y has been adjusted to provide 
the best data collapse. The same procedures have been 
used very recently in studies of unfrustrated random- 
bond Potts models lEll . 
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FIG. 6. Averaged correlation functions (normalised by 
their zero-field counterparts) against ho with y — 0.875. 
Each point is the central estimate on a strip N — 10^ sites 
long. See text and Table I for a discussion of estimated error 
bars. 



Note the use of averaged correlation functions, 
(G). We also performed plots with typical ones 
exp((ln G)), with entirely similar results. As remarked 
in Section y, one has (G) > G(ft.o — 0), on account of 
the long forward tail of the distribution. This happens 
for exp((ln G)) as well, and is a scenario valid only for low 
field intensities. Near the end of the region where scaling 
holds, on the right of Figure |^, one indeed sees the be- 
ginning of a trend towards stabilization (which would, for 
higher fields, presumably turn into a decreasing function 
of ho , were scaling still valid) . 

The value y = 0.875 used in Figure ^ gave the best 
data collapse, which still kept reasonably good over the 
interval (0.85, 0.90). The plots using exp((lnG)) behaved 
in the same way. Thus our estimate is: y — 0.875±0.025 , 
in very good agreement with the finite-size scaling ansatz 
described above, with 7 = 7/4, ly = 1. 

TABLE I. Correlation functions calculated at T — Tco, 
distance R — L and random-field intensity ho = (Go) and 
ho = 0.8L~^^^ ((G(/io))). Error bars in parentheses give un- 
certainties in last quoted digits, from spread among central 
estimates for five different runs on strips with A'^ — 10^. 
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Each point in Figure g represents an average taken 
from one run on strips N = 10^ columns long. We now 
discuss the estimation of error bars, not shown in the Fig- 
ure. Recalling that the width of the distributions is not 
expected to vanish in the thermodynamic limit, wc fol- 
low the lines extensively elaborated elsewhere for similar 
cases 1^,^^, and estimate fluctuations by evaluating 
the spread among overall averages (i.e. central estimates) 
from different samples. For values of L and ho such that 
L^/^hf) = 0.8 (approximately midway along the horizon- 
tal axis of Figure ph, we performed series of five runs, 
each with = 10 , for each L. Table I shows the re- 
sults. One sees that Equation (^ is satisfied to within 
2 parts in 10'^. Such an agreement is further evidence 
in suport of the scaling ansatz proposed above; it also 
suggests that the scaling power is y = 7/8 exactly. 

Incidentally, note that from the constancy against L of 
the ratio {G{R,L,t = 0, ho)) /G{R, L,t = 0,0), as verified 
in Table I, and the scaling of correlation functions given 
in Equation (||), one immediately has rj — riising = 1/4 
for the decay of ferromagnetic correlations at T = TcOi 
ho^O. 

We now return to scaling of the RMS relative width 
W of the distribution against field and strip width, re- 
stricting ourselves to T near Tco and ho not very large. 
For fixed R/L, taking into account that the distribu- 
tion broadens (a) with increasing random-field intensity 
(which is elementarily expected), and (b) also with in- 
creasing strip width (which we noticed in our numerics 
at r = Tco), we propose the following scaling form: 

W = h^f{Lh^) , (10) 

where the effective length Lh = hg "' plays the role of a 
saturation distance, such that f{x) constant, x ^ 1. 
In other words, (i) for high temperatures such asT — Tco 
and small ho there must be a regime in which the dis- 
tribution remains recognizably similar to Figure |^, with 
the field-induced broadening reaching a relatively small 
maximum as R,L ^ h^" (at fixed R/L). At the other 
end X <C 1, the only obvious constraint is that (ii) f{x) 
must not increase faster than x~'^^" as a; — > 0, if it does 
diverge at all. 

From scaling plots of W h^'^ against L (at T — Tco 
and ho not very large) with tentative values of the expo- 
nents, we have found the best data collapse to occur for 
K ~ 0.43 — 0.50 and u ~ 0.8. Figure 0, where the ver- 
tical axis is logarithmic, shows our results for k = 0.45 
and u = 0.8. For a; > 1, the fitting sphne is the func- 
tion y — —0.3 — 5.3 exp(— 1.57x), implying a limiting 
scaled width W h^'^ = exp(— 0.3) = 0.83 , consistent 
with (i) above. For a; < 1 the fitting curve is given by 
y = 1.73 lux — 1.40, in agreement with requisite (ii). 

To our knowledge there is no structural relationship 
between the width exponents n and u, and the stan- 
dard critical indices, such as the crossover exponent cf) 
discussed above. Conversely, one would expect widths to 
behave similarly to the above picture even at T ^ Tco, 



provided that one keeps to high temperatures and low 
field intensities. Most likely, asymptotic scaled widths 
will depend on T; a matter for further investigation is 
whether or not the numerical values of the exponents 
will also vary. 
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FIG. 7. Semi-logarithmic scaling plot of RMS relative 
widths, W ho against L . Key to symbols is the same as 
in Figure 0. Curves are fitting splines (see text). 



V. CONCLUSIONS 

We have studied the probability distributions of occur- 
rence of spin-spin correlation functions G in the d = 2 
RFIM, for binary distributions of the local fields, at 
generic distance R, temperature T and field intensity ho, 
on long strips of width L = 3 — 15 sites. 

We have shown that for moderately high temperatures, 
of the order of the zero-field transition point Tco, and 
field intensities ho ^ 0.1 — 0.2 in units of the nearest- 
neighbour coupling (the same order of magnitude used 
in most experiments), the distributions retain a rec- 
ognizable single-peaked structure, with a well-defined 
width. However, they display considerable asymmetry, 
with a short tail below the maximum and a long one 
above it, the latter owing to the mutual reinforcement 
between ferromagnetic spin-spin interactions and large 
accumulated-field fiuctuations. For low temperatures the 
single-peaked shape deteriorates markedly, as crossover 
takes place towards the double-J structure characteristic 
of the ground state. 

We have established a connection between the proba- 
bility distribution for correlation functions and the un- 
derlying distribution of accumulated field fluctuations. 
Starting from a zero-temperature description based on 
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the distribution of (essentially flat) domain walls across 
the strip, we have shown how (low-) temperature effects 
can be incorporated, and proposed analytical expressions 
for the main dependence of the distribution of correlation 
functions on R, L, T and h{). In their assumed domain 
of validity, i.e. R/L > 1, T <C 1, not very small /iq, 
and close to the upper extreme G = 1, they are in good 
quantitative agreement with numerically calculated dis- 
tributions. 

At r = Tco, for /iQ 0, we have made contact with 
scaling theory for bulk systems, and developed a finite- 
size ansatz to describe the scaling behaviour of averaged 
correlation functions. The variable that describes such 
behaviour was found to be L^Hq, with y = 0.875 ± 0.025 
from numerical data, in excellent agreement with the 
ansatz's prediction, y = 7/8. In the same region, we 
have also studied the RMS relative width W of the prob- 
ability distributions, and found that, for fixed R/L = 1 
it varies as W ^ f{L h"^) with n ~ 0.45, it ~ 0.8 . We 
have shown that f{x) fits well to a saturating form when 
X ^ oo, thus implying W ^ h!^ in d = 2. 

Further developments of the present work would in- 
clude: (i) establishing analytical expressions to connect 
field fluctuations, domain size distribution and correla- 
tion function distributions in regimes such as ~ 1 
(relevant to d = 2 behaviour), T ~ Tco, and valid for 
generic G; and (ii) a systematic study of the variation 
of widths and their associated exponents, both against 
temperature and the ratio R/L. We are currently con- 
sidering such extensions. 

Finally, as regards contact with experiment, one may 
ask how the present results for correlation functions re- 
late e.g. to the wavevector-dependent scattering ampli- 
tudes in neutron scattering [p^ . Attempts in this di- 
rection have been made earlier Since the scattering 
function reflects spatial averages over relatively extended 
regions, a connection to correlation functions must be es- 
tablished via a correlation length which rep resents the av- 
erage decay of spin-spin correlations |s|Jl8t]. Furthermore, 
fitting numerical data from one end to experimental re- 
sults from the other is a tricky task, which is usually me- 
diated by resorting to heuristically proposed line shapes. 
Of these, Lorentzian and Lorentzian-squared functions 
have been among the most popular [p|p^ , though in prin- 
ciple there is no reason why one must be restricted to 
them. A broad range of possible line shapes, compounded 
with the wide variation exhibited by several properties 
of correlation functions, as shown in the present work, 
causes one to anticipate a fairly involved investigation. 
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